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1. Introduction 



The Br avais- Pearson correlation coefficient is arguably the most widely used measure 
of dependence between random variables. For financial time series, correlations among 
returns are for instance widely used in risk management. However, there is compelling 
empirical evidence that the correlation structure of financial returns cannot be assumed 



to be constant over time, see e.g. Krishan et al. (2009). In particular, in periods of 



financial crisis, correlations often increase, a phenomenon which is sometimes referred to 
as "Diversification Meltdown". As most often potential change points are not known a 
priori, practitioners are interested in testing correlation constancy in financial time series 
at an unknown point in time. 



Wied et al. (2012) propose a nonparametric retrospective kernel-based correlation con- 



stancy test (referred to as KB-test in what follows) and Wied and Galeano (2012) propose 
a sequential monitoring procedure, complementing to other approaches for related mea- 



sures of dependence, e.g. for the whole covariance matrix (Aue et al. , 2009, Galeano and 



Pena, 2007), the copula (Na et al. 2012, Kramer and van Kampen, 2011), Spearman's 



rho ( jGaifiler and Schmid 2010), Kendall's tau (Dehling et al. 2012), autocovariances in 



a linear process (Lee et al. 2003) and covariance operators in the context of functional 



data analysis (Fremdt et al. 2012). 



In what follows, we stick to correlation. Wied et al. (2012 ) discuss why testing for constant 



correlation compared with testing for constant covariance may be of independent interest 
in applications. They also show that a correlation test can be more powerful than a 
covariance test when we have more than one change point in the covariance structure. 
However, the KB-test only considers bivariate correlations, whereas in portfolio manage- 
ment, where we typically have more than two assets, constancy of the whole correlation 
matrix is of interest. In this context, it would basically be possible to perform several 
pairwise tests and to use a level correction like Bonferroni-Holm. But probably, it is 
more elegant to directly consider the correlation matrix. This is done in the present pa- 
per which extends the methodology from the KB-test to higher dimensions, while keeping 



its nonparametric and model-free approach. We consider the 



p(p-i) 



vector of successively 



calculated pairwise correlation coefficients, derive its limiting distribution with the func- 



tional delta method approach and some proof ideas from Wied et al. (2012) and use a 
bootstrap approximation for a normalizing constant in order to approximate the asymp- 
totic limit distribution of the test statistic. This may be an alternative for the bivariate 
case as well. 



2. The fluctuation test 

Let Xf = (Xi t t,X2 t t, ■ ■ ■ ,X P)t ), t G Z, be a sequence of p-variate random vectors with 
finite 4-th moments and (unconditional) correlation matrix R t = {pt J )i<i.j< P , where 

ij _ Cov(Xi !t , Xjj) 
Pt 



vVar(X M )Var(^, t ) 

We want to test the null hypothesis H : R\ = . . . = Rt vs. H 1 : -> H . Under H , we 
denote p] 3 =: p %3 . Furthermore, we call || • || r the L r -norm, r > and D([0, 1],R ) the 
space of d-dimensional cadlag functions on the interval [0, 1], d G N, equipped with the 
supremum norm. 
The test statistic is 



V ^ k i *4j ~ij I ™ 



Q T := max \ —= \pP - pi =: max — — H-Pa-tII-, , 



where 



Pk 



J2t=i(Xi,t - Xi tk )(X jtt - Xj ;k ) 



X i,k = lJ2t=i X i,t, X J,k = l k J2t=i X 3,t an d Pk,T = (Pk ~ PT) 1<il<vUi e R P a* . The 



value p k is the empirical pairwise correlation coefficient for random variable Xi and X 



from the first k observations. Thus, the test statistic compares the pairwise successively 
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calculated correlation coefficients with the corresponding correlation coefficients calcu- 
lated from the whole sample. The null hypothesis is rejected whenever Qt becomes too 
large, that means, whenever at least one of these differences become too large over time 
or, equivalently, whenever the successively calculated correlation coefficients of at least 
one pair fluctuate too much over time. The weighting factor makes up for the fact 
that the correlations are typically estimated better in the middle or in the end of sample 
compared to the beginning of the sample. 

For deriving the limiting null distribution and local power results, some additional as- 
sumptions are necessary. The next assumptions concern moments and serial dependence 



of the random variables and basically correspond to (Al), (A2) and (A3) in Wied et al. 



(2012), adjusted for the multivariate case. 



Assumption 1. For 



( 

A i,t 



X p> t 
XitXst 



i.t) 



Hx P , t ) 

E{X ljt X 2>t ) 
E(Ai, 4 X 3 , t ) 



\Xp-ijXpj — E(Xp-i tt X p>t ) J 



and Sj := Yjt=i Ut, we have 



lim E ( -S m S' m ) =: D lt 



where D% is a finite and positive definite matrix with 2p + 



p(p-i) 



rows and 2p + 



p(p-i) 



4 



columns. 

Assumption 2. For some r > 2, the r-th absolute moments of the components ofU t are 
uniformly bounded, that means, sup t6Z E||[/ t || r < oo. 



Assumption 3. The vector . . . , X Pjt ) is L 2 -NED (near-epoch dependent) with size 
— r ^z^ with r from Assumption^ and constants (c t ),t G Z, on a sequence (V t ),t G Z ; 
which is a-mixing of size 4>* : = — ^2, i-e. 

||(X M , . . .,X Ptt ) - E ((A ljt , . . .,X p , t )\<r(Vt-i, ■ V t+l ))\\ 2 < c tVl 

with lim^oo Vi = such that c t < 2||f/ t || 2 with U t from Assumption^ 

Assumption [T] is a regularity condition which holds in most models, except perhaps the 
case of trending random variables. This situation is however not relevant for financial 
returns. 

Assumption [2] is more critical because it requires finite |4 + <5|-th moments of X t with 
5 > arbitrary (note that the components of X t enter U t quadratically). In fact, there 



is evidence that even variances might not exist for some financial series, cf. Mandel 



brot (1962). However, simulation evidence shows that the test still works under the 
t3-distribution (see below). 

Assumption |3] is a very general serial dependence assumption which holds in most relevant 



econometric models, e.g. in GARCH-models under certain conditions (cf. Carrasco and 



Chen, 2002). It guarantees that the vector 

0^i,t> • • • > ^p,ti -^i,t> • • • 1 X P ,t, X\^2,ti XifX^j, . . . , Xp-ijXpj 



is L2-NED (near-epoch dependent) with size — |, cf. Davidson (1994), p. 273, which is 



necessary for applying a functional central limit theorem later on. 



Next, we impose a stationarity condition which is in line with Aue et al. (2009). 



Assumption 4. Under Hq, (X^t, ■ ■ ■ ,X Pit ),t G Z, is weak-sense stationary. 



This condition might be slightly relaxed to allow for some fluctuations in the first and 



second moments (see (A4) and (A5) in Wied et al. 2012), but for ease of exposition and 
because the procedure would remain exactly the same, we stick to this notation. Note 
that most financial time series processes as for example GARCH are (unconditionally) 
stationary under certain conditions. Clearly, the original test problem is invariant under 
heteroscedasticity, but we believe that it is at least very much involved if not impossible to 
design a fluctuation test for correlations in which arbitrary variance changes are allowed 
under the null hypothesis. 
Our main result is: 



Theorem 1. Under Hq and Assumptions l^^J 



(Pt( s ) ~ Pr)i<i,3<P,irt =^ E l/2 B 



p(p-i) 



P(P-1) w P(P~1) 

2 X 2 



on D (J0,1],M^) ; where t(s) = [2 + s(T-2)] ; E ~ ( J - 
depends on the data generating process and B v(v ^ ' (s) is a vector of p ( p ~^ independent 
standard Brownian bridges. 

The proof of the theorem is in the appendix. It relies on the application of an adapted 
functional delta method. We want to stress that simply applying a functional central 
limit theorem is not enough here due to the cumbersome, non-linear structure of the 
correlation coefficient. 

To obtain critical values, we need information about E. There are several possibilities for 
estimating E; one possibility is the estimator E, given by the bootstrapped covariance 
matrix of the vector \/T (pr)i<i j< P i^j- ^ or ^ s estimation, one can for example use the 



moving block bootstrap from Kiinsch (1989) and Liu and Singh (1992), cf. also Lahiri 



(1999), Congalves and White (2002), Calhoun (2011), Radulovic (2012) and Sharipov 



and Wendler (2012). 



Defining a block length It, we divide the time series into T — — 1 overlapping blocks 



Bi, i = 1, . . . , T— l t — 1, with length It such that B\ = (Xi, . . . , X/ T ), B 2 = (X 2 , . . . , X; T+ i), 



times 



Then, in each bootstrap repetition b,b = 1, . . . , B for some large B, we sample 
with replacement one of the T — It — 1 blocks and stick the blocks together so that we 
obtain B p-dimensional time series with length ^ • It- For each bootstrapped time 
series we calculate the vector VT (pJrp) . . .. The estimator E is then the boot- 
strap covariance matrix of these B vectors, that means the variance with respect to the 
bootstrap sample conditional on the original data. 

The bootstrap estimator "replaces" the rather complicated kernel estimator E from the 



KB-test (Appendix A.l in Wied et al. , 2012). The advantage of the bootstrap estimator 



is the fact that it can be derived easily even in higher dimensions. It would basically be 
possible to obtain a kernel estimator also in higher (> 2) dimensions, but its structure 
would then depend on the structure of derivatives of certain higher-dimensional func- 
tions (cf. the proofs) which makes it very cumbersome and much harder to implement. 
Moreover, a kernel estimator depends on the choice of the bandwidth and the kernel. 
The disadvantage of the bootstrap is that it is computationally more intensive and that 
stronger moment conditions are required for the proofs (see below). In addition, the 
choice of the block length is required. 

For asymptotic results concerning the bootstrap, we need an assumption on the block 
length. 

Assumption 5. For T — > oo, Z r — » oo and l T = o(T). 



This assumption comes from Calhoun (2011), Corollary 2; it guarantees that the block 



length becomes large but not too large compared to T. 



Theorem 2. Under H and Assumptions T^^^^S, a consistent estimator for E is given 
by the bootstrapped variance (block bootstrap) of the vector y/T (fir)^ <pi ^ 

The following corollary yields the asymptotic null distribution of the test statistic itself. 



Corollary 1. Under Hq and Assumptions l^^l^S 



Qi 



max 

2<k<T 



Vr o< s <i 



E^B 



p(p-i) 



and 

k 



max 



E' 1 ' 2 P l 



k,T 



->d sup 

1 0<s<l 



p(p-i) , . 

B 2 (s) 



2<k<T y/T 

There might be situations in practice in which E is not positively definite so that E~ x l 2 
would not be defined. However, due to Assumption [T], at least for larger T, we can 
virtually assume positive definiteness. 

The quantiles of the limit distribution from Corollary [TJ which serves as an approximation 
for the finite sample distribution, can easily obtained by Monte Carlo simulations, that 
means by approximating the paths of the Brownian bridge on fine grids. 

3. Local power 

Statisticians are often not only interested in the behavior of a test under the null hypoth- 
esis, but wish to get information about the behavior under some local alternatives. For 
simplicity, we consider a setting in which the expectations and variances remain constant 
such that a covariance change is equal to a change in correlations. To derive specific 
results, we impose 

Assumption 6. Under Hi, there is a correlation change of order ^= (M > arbitrary) 
in at least one of the components of X 4 with constant expectations and variances and 



{E( X i,tXj,t))i<i,j< P ,i^j - k + -j=g ^— J 



where k £ IR P(P 2 1 is a constant vector and g(s) = (gi(s), . . . , gr v<.v-i) \(s)) is a bounded 
p( - p ~ 1 - > -dimensional function that is not constant and that can be approximated by step 
functions such that the function 

g{u)du — si g{u)du 
o Jo 

is different from G 1R P<P 2 for at least one s G [0,1]. Here, the integral is defined 
component by component. 



Note that we now basically deal with triangular arrays because the distribution of the 
X t changes with T, but for simplicity we do not change our notation. 
A typical example for the function g would be a step function with a jump from to 
go in a given point Zq in one of the components which implies that the correlation of 
one pair jumps at time [T • zq]. A step function with several jumps would correspond 
to multiple change points. Using a continuous function g would lead to continuously 
changing correlations. 

The following Theorem [3] is an analogue to Theorem [T] and yields the distribution under 
the sequence of local alternatives. 



Theorem 3. Under Hi and Assumptions ~T^§^6_ 



t{s) 



r(s) 



where 



C(s) =M 



" Q s gi{u)du - s $1 gi{u)du 



, 1 ( q 2 (u)du — s L q 2 (u)du 



i / 1 ( fn Qp(p-i) (u)du — s L Op(p-i) (u)du ) , 

is a deterministic function that depends on the specific form of the local alternative under 
consideration, characterized by g. 

In Theorem [3j the supremum is taken over the absolute value of a Brownian bridge plus a 
deterministic function C(s). The main characteristic of the function C(s) from Theorem 
[3] is that we have in each component i = 1, . . . , p< - p ~ 1 ^ ; the factor M times the expression 
f Q S 9i(u)du — s Jq 1 gi(u)du, which follows from the structure of a Brownian bridge. 
Also under local alternatives we want to estimate E with the bootstrap and it turns 
out that the estimator presented in Section [2] has the same limit distribution as under 
the null hypothesis. Thus, the bootstrap approach (which does not require weak-sense 
stationarity) is valid both under the null and under the alternative. 



Theorem 4. Under Hi and Assumptions 7|£|3][5][6, a consistent estimator for E is given 
by the bootstrapped variance (block bootstrap) of the vector y/T {Pr)i<i < P i ^ ;■ 

Corollary 2. Under Hi and Assumptions 7]|£|[3|5][6, 



Qt = max —=\\P kT \\, -» d sup 

2<fc<T y/T 0<s<l 



E l l 2 B Pj£ ^{s) + E l l 2 C{s) 



and 



k 



max 



E~ l ' 2 P 



k,T 



->-d sup 

1 0<s<l 



p(p-i) 



2<fc<T ^/J 1 

Corollary [2] basically serves for two aspects. First, for known g (and thus C) we get 
information about the distribution of the limit under the local alternative. Moreover, the 
result serves for showing "consistency" of the multivariate correlation test as the following 
corollary shows. 

Corollary 3. Under assumptions l^^S^B, the fluctuation test is consistent in the follow- 
ing sense: If Ph 1 {M) denotes the rejection probability for given M under the alternative 
and e > is fixed, then there is an Mq such that lim^oo Phx{M) > 1 — e for all M > Mq. 



4. Finite sample evidence 

We illustrate the finite sample properties of our multivariate test with Monte Carlo 
simulations in different settings: We consider a series of four-variate random vectors 
which are, on the one hand, serially independent and, on the other hand, fulfill a four- 
variate MA(l)-structure with MA-parameters 0.1. The lengths of the series are chosen 
as T = 200,500,800,1000, the block lengths are It = [T 1 ^ 4 ], respectively^] the number 
of bootstrap replications is 199 and the number of Monte Carlo replications is 1000. We 
consider, on the one hand, a four-variate normal distribution (ND) and, on the other 
hand, a four-variate ^-distribution. The ^-distribution is not covered by our assump- 
tions, but we analyze it to get a picture of the behavior of the test in settings which are 



iFor, T = 200, 



j- , so that the length of the bootstrapped time series is not exactly equal to 



T . However, we consider the difference as negligible. 
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realistic in financial applications. For simulating the behavior under the null, we consider 
constant variances 1 and pn = . . . = pu = =: p. For the behavior under alternatives, 
we consider a change in the correlation between X\ and X2, on the one hand, and a 
change in all six pairwise correlations, on the other hand, in the middle of the sample 
with shifts Ap = -0.2, -0.4, -0.6, 0.2, 0.4, 0.6. 

The results (empirical rejection probabilities, not-size-adjusted, nominal level 0.05 which 
corresponds to a simulated critical value of 4.47) are given in Table [XJ 

— Table Q] here - 

It is seen that the size is basically kept, even for the heavy-tailed distribution. The power 
of the test increases in T and in absolute values of the correlation changes and is much 
higher if all correlations change. Interestingly, a decrease of all correlations is better 
detected than an increase of all correlations; this is not true in case of a change in only 
one component. For the t 3 -distribution, the power is in general considerably lower. For 
serial dependence, the results basically remain the same. 

In another setting, we compare the bivariate bootstrap with the bivariate kernel-based 
test. The setting is basically similar to the first setting. We choose the Bartlett kernel 
and [log(T)\ as bandwidth for the kernel-based test and consider the cases T = 200, 500. 
The results (empirical rejection probabilities, not-size-adjusted, nominal level 0.05 which 
corresponds to a critical value of 1.358) are given in Table [2j 

- Table [2] here - 

We see that in general, the results are more or less the same; the small differences may 
partly be caused by Monte Carlo inaccuracies. We repeated these simulations also for 
heavier serial correlation and as expected we had strong size distortions in both cases. 
However, these distortions can be avoided by using a larger bandwidth or a larger block 
length, respectively. 

Furthermore, we observe that the power of the bivariate test for a certain change Ap is 
larger than the power of the four-variate test when there is a change in just one pair, but 
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lower than the power of the four-variate test when there is a change in all six pairs. This 
is plausible due do the construction of the multivariate test. 



5. Application to stock returns 

Next, we show how the proposed test can be applied in financial time series by considering 
the correlation of four stocks. In order to avoid issues due to market trading in different 
time zones, we just consider the European market. We look at the four companies of 
Euro Stoxx 50 with the highest weights in the index in the end of May 2012, that means 
Total, Sanofi, Siemens and BASF, and consider the time span 01.01.2007 - 01.06.2012 
such that T = 1414. The data was obtained from the database Datastream. Figures [TJ 
[2j [3] plot rolling windows of the six pairwise correlations of the continuous daily returns 
from each asset with the window length 120, corresponding to the trading time of about 
half a year. The days on the x-axis show the first day of the windows, respectively. 

- Figure [l] here - 

- Figure [2] here - 

- Figure [3] here - 

We identify time- varying correlations; it is for example interesting to see that the corre- 
lation between Total and Sanofi is close to in the beginning of February 2008 and much 
higher after this. The correlation between Sanofi and BASF is interestingly low in the 
middle of 2009. 

The test statistic Qt applied on the four-variate return vector is equal to 10.49. With this 
value, we cannot yet determine if the null hypothesis of constant correlation is rejected. 
So we calculate the term max 2 <KT -ym E~ l ^ 2 Pk t , which has to be compared with 

VT > x 

the quantiles of sup 0<s<1 ||_B 6 (t)|| 1 , see Corollary [TJ with different amounts of bootstrap 
replications (but the same starting value, respectively). The results are given in Table [3j 
We see that the values slightly fluctuate, but, as expected from the theory, converge to a 
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proper limit which yields the result that the null hypothesis is rejected on the significance 
level a = 0.05. The approximate p-value is smaller than 0.001. 
Figure [4] shows the process 

\l<i,3<P,i^3 V / 2 <k<T 

that means the evolution of the successively calculated correlations over time. In the 
context of CUSUM tests, the point of the maximum is often considered as a reasonable 
estimator for the (most important) change point if the test decides that such a point 



actually exists, see Vostrikova (1981) and the related literature. In this case, the max- 
imum is obtained in the middle of September 2008, corresponding quite closely to the 
insolvency of Lehman Brothers. 

- Figure [4] here - 



6. Conclusion 

We have presented a new fluctuation test for constant correlations in the multivariate 
setting for which the location of potential change points need not be specified a priori. The 
new test bases on a bootstrap approximation, works under mild assumptions regarding 
the dependence structure and has appealing properties in simulations and on empirical 
data. It might be an interesting question for the future to extend the present approach to 
the problem of monitoring correlation changes or to other, perhaps more robust measures 
of dependence. 
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A. Appendix 



Proof of Theorem [7] 

At first, we need an invariance principle for the vector 

X p,t - 

X Pi t — E(X P) t) 

X\ t tX 2j t — E(Xi )t X 2) t j/ 

X\f,X?,£ — E^XijXzj) 

\X p -ijX p ^ — E(Xp„i i tX P)t ) i 



1 r(s) 

v A t=i 



E(4 



E(X lit ) 



which is provided by 



Davidson 



on D([0,1],R 2 p+^Y where W^+^fV) is a f *&=^ 



jl994[), p. 492. Thus, it holds V T (s) ^ d D\ /2 W 3p+sSs r !1 {s) 

+ 2p j- dimensional Brownian 



motion with independent components. 
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Now, noting that 



V T (s) :-- 



( X\ 



X p 



X1X2 — 



X1X3 



E(X M ) 

E(X Rt ) 
E(X lit X 2)t ) 
E(Xi,iX 3 , t ) 



\Xp-iXp — E(Xp_i jt X P)f ) y 

the goal is to transform this vector of simple first and second order moments into the 
vector with the successively calculated correlation coefficients and then to apply the 



adapted functional delta method, Theorem A.l in Wied et al. (2012). The transforming 
functions are 







p(p-l) 

2 






( 


x x 


( X p+l) 








( X 2p) 


(^•••^( 2p+ Mp-i))) -> 




x 2p+l 


XpXp+i 






X2p+2 


XpXp+2 




V( 


2 P+ ^- 1) ) 


— X2p-\X2p j 
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for the transformation on the vector of variances and covariances and 



. „ i p(p-i) _ p(p-i) 



/ Zp+1 \ 



(ari, . . .,x/ p| p(P -i) Q -)■ 



for the transformation on the vector of correlations. 
We obtain 



W T (s) :-- 



(p 



t(s) 



p J h< 



p(p-i) , . 



j2p+ p(p 2 ' j £ Qr m atrices Do 



p + &£±) x (2p + 



p(p+i) 

2 



(1) 



and Dz 



on £>([e,l], 

^ p(p+i) x ^ _|_ p(p+i) j j Here, L*2 is the Jacobian matrix of /1 and -D3 is the Jacobian 
matrix of f 2 , evaluated at certain moments. For arbitrary e > 0, D Me, 1] 



p(p-i) 



the space of (2p + p( - p ~ 1 - > j -dimensional cadlag functions on the interval [e, 1], equipped 
with the supremum norm. 

We are not interested in the exact (and cumbersome) structure of these matrices, but we 
observe that D 2 contains all [p + p ( p ^ j -dimensional unit vectors and D3 contains all 
^ g(p±jO j -dimensional unit vectors (weighted with some constants) in its columns. Thus, 
D 2 and D 3 have full column rank. 



Now, with an application of Theorem 4.2 in Billingsley (1968), we obtain Wt(s) =^d 
DzD 2 D\ /2 W 2ip+S ^ (s) on D ( [0, 1], R^+^^Y Moreover, it holds 



D 3 D 2 Dl /2 W 2p+Pj ^(s) 1 (DsDzD^D'^W^is) 



and then, the proof follows with the continuous mapping theorem. 
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Proof of Theorem 

From it is easy to see (with s = 1) that the asymptotic covariance matrix of 
^ (/'r),. , r is ec l ual t0 D 3 D 2 D 1 D' 2 D' 3 =: E. 



Now note that the bootstrap theorem for near epoch dependent data from Calhoun (2011 ), 
Corollary 2 is applicable on £/r(l) using the standard Cramer-Wold device for the mul- 
tivariate generalization. Assumption 1 is fulfilled with our Assumption |3j Assumption 2 
with our Assumption [4], Assumption 3 with our Assumption [2] and Assumption 4 with 
our Assumption [5j So, we can consistently estimate the asymptotic covariance matrix of 
{7r(l) with the block bootstrap. But then, with the standard (functional) delta method 



for the bootstrap (van der Vaart and Wellner, 1996 Theorem 3.9.11.) transforming C/y(l 



to Vt(1) the validity of the bootstrap directly follows. ■ 
Proof of Corollary [7] 

This follows with the continuous mapping theorem noting that D is positively definite 
due to Assumption [T] and given the remarks regarding the regularity of D 2 and D% from 
the proof of Theorem [TJ ■ 

Proof of Theorem [3] 

Transferring the proof of Theorem[lJ we obtain, under Hi, Vy(s) =^ D^ 2 W 2p+P{P 2 " (s) + 
A(s) on D (jO,l],M 2p+ ^), where 



(note that g(u)' is the transpose of the function g). 
So, 



Wt(s) := ^(P% } ~ //•'),.,,■ „,•„• ^ D^W PJ ^(s) + D 3 D 2 A(s), (2) 



where D 3 and D 2 are the matrices mentioned in the proof of Theorem [T] Due to the 
structure of D 3 and D 2 , we have 
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D 3 D 2 A(s) = M 



-, 1 f n s q 1 (u)du 

v /Var(Xi)Var(X 2 ) J ° ' 

1 ff q 2 (u)du 

v /Var(X 1 )Var(X 3 ) J ° WA ' 



\ 



i / 1 fn 9p(p-1) (u)du , 



This completes the proof. 



Proof of Theorem 



The local alternatives fulfill Assumption 2 in Calhoun (2011), so that, by Corollary 2 in 



Congalves and White (2002) and the Cramer- Wold device, we estimate E consistently 



with the bootstrap estimator. 
Proof of Corollary [£] 

This follows in the same way as in the proof of Corollary [T} 



Proof of Corollary [3] 

For large M the test statistic becomes arbitrarily large, especially larger than every quan- 
tile of the limit distribution under the null hypothesis. ■ 
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Table 1: Empirical size and empirical power (times 100, respectively) of the multivariate corre- 
lation test; column 3-7 gives empirical rejection probabilities for a change in correlation between 
X\ and X% column 8-9 gives gives empirical rejection probabilities for a change in correlation 
between X\,X2-,X% and X4, respectively (there are no results for * because in these cases, the 
correlation matrices are not positively definite) 



MA 


distr. 


Ap 


rej.prob. 


rej.prob. 








T = 200 


T = 500 


T = 800 


T = 1000 


T = 200 


T = 500 





N 





2.7 


5.5 


5.6 


6.4 


4.3 


5.3 





N 


0.2 


8.1 


16.3 


30.6 


33.1 


36.5 


87.7 





N 


0.4 


21.4 


56.8 


79.1 


92.9 


93.1 


100.0 





N 


0.6 


38.5 


87.2 


98.5 


100.0 


100.0 


100.0 





N 


-0.2 


7.0 


18.4 


28.4 


35.0 


77.6 


100.0 





N 


-0.4 


22.8 


56.6 


81.8 


91.7 


* 


* 





N 


-0.6 


37.1 


86.8 


98.9 


100.0 


* 


* 





t 





12.1 


8.5 


6.2 


6.8 


10.0 


8.9 





t 


0.2 


13.5 


13.1 


14.1 


14.9 


20.1 


35.2 





t 


0.4 


20.4 


24.9 


33.4 


38.9 


47.8 


81.2 





t 


0.6 


28.9 


45.8 


56.0 


65.9 


77.3 


95.3 





t 


-0.2 


13.2 


13.7 


13.7 


14.5 


40.0 


66.0 





t 


-0.4 


18.7 


26.0 


34.4 


38.3 


* 


* 





t 


-0.6 


31.8 


49.5 


57.9 


65.2 


* 


* 


0.1 


N 





4.1 


5.9 


6.3 


5.2 


3.9 


4.6 


0.1 


N 


0.2 


7.8 


19.5 


29.5 


33.1 


37.3 


89.8 


0.1 


N 


0.4 


20.1 


55.1 


78.2 


88.2 


92.0 


100.0 


0.1 


N 


0.6 


40.1 


87.4 


98.8 


99.7 


100.0 


100.0 


0.1 


N 


-0.2 


8.8 


18.1 


27.3 


35.7 


74.8 


100.0 


0.1 


N 


-0.4 


21.5 


55.6 


77.9 


87.3 


* 


* 


0.1 


N 


-0.6 


38.5 


91.6 


99.2 


100.0 


* 


* 


0.1 


t 





11.1 


8.7 


6.8 


6.8 


11.6 


9.7 


0.1 


t 


0.2 


11.5 


13.4 


13.0 


14.6 


21.5 


36.4 


0.1 


t 


0.4 


20.3 


25.3 


33.1 


36.8 


49.3 


83.1 


0.1 


t 


0.6 


29.4 


45.8 


59.0 


64.9 


77.7 


96.1 


0.1 


t 


-0.2 


13.1 


14.1 


13.3 


16.4 


37.7 


63.9 


0.1 


t 


-0.4 


20.3 


24.5 


31.2 


40.2 


* 


* 


0.1 


t 


-0.6 


29.7 


44.0 


60.2 


65.3 


* 


* 
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Table 2: Empirical size and empirical power (times 100, respectively) of the bivariate correlation 
test; KB denotes the kernel-based test and BS denotes the bootstrap test 



MA 


distr. 


Ap 


rej.prob. 








T = 200, KB 


T = 200, BS 


T = 500, KB 


T = 500, BS 





N 





4.1 


4.4 


3.1 


6.1 





N 


0.2 


20.9 


21.8 


47.2 


49.7 





N 


0.4 


66.1 


70.4 


98.3 


97.8 





N 


0.6 


96.9 


96.7 


100.0 


100.0 





N 


-0.2 


19.5 


21.1 


48.3 


45.3 





N 


-0.4 


65.2 


69.6 


98.6 


97.3 





N 


-0.6 


97.2 


96.4 


100.0 


100.0 





t 





7.0 


6.1 


4.7 


3.7 





t 


0.2 


10.6 


10.5 


13.1 


13.3 





t 


0.4 


26.4 


29.8 


43.5 


46.0 





t 


0.6 


52.2 


54.6 


74.8 


76.6 





t 


-0.2 


11.7 


11.4 


14.0 


12.9 





t 


-0.4 


28.1 


28.1 


41.9 


44.6 





t 


-0.6 


52.1 


53.7 


77.0 


77.7 


0.1 


N 





4.0 


5.1 


3.6 


4.2 


0.1 


N 


0.2 


20.4 


19.5 


47.9 


47.9 


0.1 


N 


0.4 


68.1 


68.3 


98.0 


98.8 


0.1 


N 


0.6 


97.0 


97.1 


100.0 


100.0 


0.1 


N 


-0.2 


19.5 


20.5 


49.6 


45.7 


0.1 


N 


-0.4 


68.9 


68.1 


98.6 


98.9 


0.1 


N 


-0.6 


96.5 


96.9 


100.0 


100.0 


0.1 


t 





6.2 


5.5 


3.9 


3.8 


0.1 


t 


0.2 


11.0 


11.1 


13.1 


15.8 


0.1 


t 


0.4 


27.4 


27.4 


46.0 


46.5 


0.1 


t 


0.6 


54.2 


55.5 


72.8 


77.2 


0.1 


t 


-0.2 


11.2 


11.0 


13.3 


15.3 


0.1 


t 


-0.4 


28.2 


28.2 


45.5 


44.3 


0.1 


t 


-0.6 


52.6 


55.0 


72.4 


74.7 



Table 3: Standardized test statistic in the empirical illustration 

b 199 599 1099 2099 5099 10999 critical value, a = 0.05 

D- 1/2 Q T 6.71 6.37 6^36 646 6^51 6^55 4.47 
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2007 



2008 



2009 2010 
Year 



2011 



(a) Rolling correlations between Total and Sanofi 




2007 



2008 



2009 2010 
Year 



2011 



(b) Rolling correlations between Total and Siemens 
Figure 1: Rolling correlations 
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2007 2008 2009 2010 2011 

Year 

(a) Rolling correlations between Total and BASF 




2007 2008 2009 2010 2011 

Year 



(b) Rolling correlations between Sanofi and Siemens 
Figure 2: Rolling correlations 
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hi 

2007 2008 2009 2010 2oTT 

Year 

(a) Rolling correlations between Sanofi and BASF 




2007 2008 2009 2010 2011 

Year 



(b) Rolling correlations between Siemens and BASF 
Figure 3: Rolling correlations 
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2007 2008 2009 2010 2011 2012 

Year 

Figure 4: Evolution of successively calculated correlations 



27 



